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ABSTRACT 

The rotation rates in the deep interior and at the surface of 22 main-sequence stars with masses 
between 1.0 and 1.6 Mq are constrained by combining asteroseismological analysis with 
spectroscopic measurements. The asteroseismic data of each star are taken by the Kepler or 
CoRoT space mission. It is found that the difference between the surface rotation rate and 
the average rotation rate (excluding the convective core) of most of stars is small enough to 
suggest that an efficient process of angular momentum transport operates during and/or before 
the main-sequence stage of stars. If each of the surface convective zone and the underlying ra¬ 
diative zone, for individual stars, is assumed to rotate uniformly, the difference in the rotation 
rate between the two zones turns out to be no more than a factor of two in most of the stars 
independently of their ages. 

Key words; asteroseismology, methods: data analysis, stars: interiors, stars: oscillations, 
stars: rotation, stars: solar-type 


1 INTRODUCTION AND SUMMARY 


Rotation is one of the fundamental issues that affect stellar struc¬ 
ture and evolution. It induces material mixing processes which may 
have a significant influence on nucleosynthesis in stars and affects 
their evolution tracks. In addition, the stellar rotation plays a crucial 
role in the dynamo process, which is essential for the generation 
and maintenance of the magnetic field in stars. 

A key process to understand the stellar rotation is angular 
momentum transport. In almost all phases of evolution, the core 
of stars gradually contracts, while the envelope expands, so that 
the rotation rate of the core is larger than that of the envelope, as 
long as the angular momentum is locally conserved. This general 
trend is compensated by the transport of the angular momentum 
from the core to the envelope. While many uncertainties still re¬ 
main in the theoretical treatment, some clear evidences to support 
the presence of efficient angular momentum transport have been 
found in observations. Helioseismology has demonstrated nearly- 
uniform rotation of the Sun in the radial directi on from 30 per cent 
of the total radius (at l east) to the surface (e.g. Scho^rt 3 19981: 
iThompson et alj|2003l : lEff-Darwich & Korzennild 20131) . The sur¬ 
face rotation rates of white dwarfs, which are essentially cores of 
red giants, are much smaller than the value obtained by shrinking 
the Sun to the typical size of those stars assuming local conser- 
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vation of angular momentum. Asteroseismology has also recently 
shown that the core-to-envelope ratio of the rotation rate of sev¬ 
eral subgiants and young red giants is about 60 at most dBeck et alj 
l2012l : lDeheuvels et al.ll2012Ll2014l) . and that almost uniform rota¬ 
tion is detected in a main-s equence A star ^ Kurtz et ^l2014h and 
in a main-sequence F star dSaio et alj|2015 ). All of these results 
indicate that the angular momentum in the core is efficiently trans¬ 
ported to the envelope. On the other hand, we remark that some 
model-dependent studies about B stars based on ground-based ob- 
servation have presented only weak arguments of non-rigid rotat ion 

dAerts et alj|2003l : IPamvatnvkh. Handler & Dziembowskill2004t) . 

In spite of these results, the number of main-sequence stars in 
which the structure of the internal rotation has been measured is 
still very limited. In this paper, we present the first statistical study 
of the radial rotation structure for main-sequence solar-like stars. A 
method is proposed to constrain the ratio between the internal and 
surface rotation rates of a star that shows solar-like oscillations. On 
one hand, we measure the average interior rotation rate, /sets, by as¬ 
teroseismology. On the other hand, the surface rotation rate, /surf, 
is determined either by spectroscopy or by measuring periodic lu¬ 
minosity variation due to stellar spots. An agreement between /seis 
and /surf implies uniform rotation, whereas a disagreement sug¬ 
gests a differential rotation. The determination of the core rota¬ 
tion rate based on high-order and low-degree p modes has been 
tried in helioseismology (e.g. lElsworth et alj 1 19951 : lUazrek et alj 
ll996l:lGarcfa~eF ^120041) . though these have been found to be only 
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weakly sensitive to the s olar rotation bellow 0.2Rq dChaplin et alj 
ll999l:[Garciaetal.l2008h . Comparison of the core rotation rate with 
the surface rotation rate derived from luminosity variation due to 
spots has been suggested in the case of solar-lik e pulsators by, 
e.g., lAerts, Christensen-Dalsgaard & Ku^ l l201(]li . In this paper, 
we consider not only spots but also spectroscopic measurements 
to constrain the surface rotation rate, and find that the latter is more 
reliable. 

The method has been applied to a sample of main-sequence 
solar-lik e pulsators, 20 of wh ich have been observed by the Kepler 
mission j Borucki et al.ll201^ . and two by the CoRoT spacecraft 
jBaglin et alJl2006iJlbh . These are 1.0-1.6 Mq stars with the solar 
composition, while their evolutionary stages cover the entire main- 
sequence phase. None of these stars are known to be members of a 
close binary system. In addition, they rotate faster than the Sun by 
(typically) an order of magnitude. 

The comparison of the surface rotation rate with the average 
internal rotation rate shows that the two rates are generally close 
to each other. For 10 stars, their difference is too small to be ex¬ 
plained by simple evolutionary models that assume local conser¬ 
vation of angular momentum in the radiative layers, and constant 
rotation rates in the convection zones. This clearly demonstrates 
that angular momentum is transported from the inner layers to the 
outer layers (in the radiative layers and/or at the interfaces between 
the radiative and convective zones) of these stars. The result would 
generally be valid for main-sequence stars with similar masses to 
the Sun. Under the assumption that the surface convective zone and 
the internal radiative zone rotate rigidly (with different rates from 
each other in general), we found that 21 out of 22 stars show rota¬ 
tion rates in the radiative interior that are consistent with those in 
the convective envelope within a factor of two. 

The paper is structured as follows, sectionj^presents how the 
radial structure of the rotation profile can be measured in solar-like 
stars. Then section [3 describes the observational constraints and 
the method, followed by section |4] that shows our results. Finally, 
sectionj^is devoted to discussion. 


2 ASTEROSEISMIC MEASUREMENT OF INTERNAL 
ROTATION 


2.1 Mode properties 

Asteroseismology is a powerful method to probe physical proper¬ 
ties of invisible stellar interiors, using pulsations detected at the sur¬ 
face. Recent space-borne observations made by the CoRoT and Ke¬ 
pler satellites have detected oscillations of a la rge number of solar¬ 
like s t ars with an unpreced ented precision (e.g Appourchaux et al] 
I2OI2I: iMathuretalJ 120121) . opening the possibility of scrutinis¬ 
ing the internal structure o f the stars similar to the Sun statis¬ 
tically iChaplin et alj|20f4h . More importantly, the rotation rates 
both at the surface and in the deep layers have been measured 
for a significant number of stars at various evolutionary stages 
(e.g. [Be nomar. Appourc hau x & Ba udinll2 TO9l: iBarban et a.l.|[20()9l : 
lOeheuvels et alj 2012L I20l4 lOavies et alj 2015t) . The measure of 
the internal rotation is made possible because rotation induces a 
splitting of oscillation eigenfrequencies, which could be compared 
to the Zeeman effect observed in spectral lines in presence of a 
magnetic field. More precisely, for a non-rotating star, a pulsation 
mode is described by two integers n and I, which correspond to 
the number of nodes along the radius of a star (radial order of the 


mo d4il) and the number of nodal lines on the surface (degree of the 
mode), respectively. There exist 2/ -f 1 modes with the same n and 
I that have a common eigenfrequency (degenerated solutions), but 
different dependence on the azimuthal angle, if the star is spheri¬ 
cally symmetric. However, the rotation lifts the degeneracy among 
non-radial pulsation modes (I > 0 ), revealing a discrete structure of 
modes identified by their azimuthal order m such as —Z ^ m ^ -\-l, 
with \m\ equal to the number of nodal lines on the surface that co¬ 
incide with meridians. Slowly rotating stars like the Sun almost 
preserve their spherical symmetry^ so that each mode frequency 
is given by 

I'n.i.m = -f m5Vri,l- (1) 

Here, Vn,i is the frequency for a spherically symmetric star, while 
5vn,i is called the rotational splitting, which represents a small per¬ 
turbation induced by the rotation. The rotational splitting is the ob¬ 
servable parameter that enables us to measure the rotation of the 
stellar interior. 


2.2 Rotational kernels 


We explain how the rotational splittings are sensitive to the internal 
rotation profile of the star, assuming for simplicity that the rotation 
rate is a function of only the distance from the centre of the star. 
The sensitivity depends on the properties of observed oscillation 
modes, which are high-order (large n) and low-degree (small Z) p 
modes consisting of acoustic waves, in the case of solar-like pul¬ 
sators. The constituent waves can travel from the surface to deep 
layers before reaching their turning point, where they are refracted 
back toward the surface. The sensitivity of such modes to the rota¬ 
tion profile within the visited layers is represented by the so-called 
rotational kernels, hereafter simply referred as kernels, Kn,i{r). 
The kernel of each mode mostly depends on the eigenfunction and 
is related to the rotational splitting, 60 ^, 1 , and to the cyclic fre¬ 
quency of rotation, /(r), by 

5tyn,i = (1-C„,i) f Kn,i{r)f{r)dr 

Jo 

= ( 2 ) 


Thus the splitting is proportional to the average cyclic rotation fre¬ 
quency fn,i, across the region visited by a mod e of degree Z and ra¬ 
dial order n. Here, C n.i is the Ledoux constant dCowling & Newing 
ll949l : lLedouxll95l]) . This constant measures the effect of the Cori¬ 
olis force and tends toward 0 for high-order p modes. For example, 
Cn,i — 10“^ in the range of observed modes in the Sun. This im¬ 
plies that the splitting may differ from the rotation frequency 
by about 1 per cent at most. Because the typical uncertainty in the 
rotational splittings of solar-like stars is as large as 5 per cent, the 
Coriolis term (7^,; is in practice negligible in equation 

Figs [Ha) and [^a) show examples of kernels for n = 10 
and n — 25 for the models of the three different stars; the 
Sun ( IChristensen-Dalsg aaMeljilj[jj96l). Kepler-25 (M — 1.26 ± 


0.03 M 0 . see lBenomar et alj|2014bl) and HAT-P-7 (M = 1.59 ± 


^ This is only an approxim ate explanation about the radial order. More 
detailed discussion is found in lXakatal j2012h . 

^ Strictly speaking, the effect of the centrifugal force that causes distortion 
of the stellai' configuration might not always be negligible in the stars dealt 
in this paper, because they rotate significantly faster than the Sun. This prob¬ 
lem is discussed in section[33] On the other hand, we assume for simplicity 
that the effect of the magnetic field is much smaller than that of rotation. 
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0.03 M 0 . iBenomar et alj2014bl:lLund et alj20l4) . Only modes of 
degrees I = 1 and I = 2 are shown, as it is rare to observe / = 3 
(and higher-degree) modes in Kepler or CoRoT data. The radius 
is normalised by the stellar radius R in each case. The kernels are 
denser and of greater amplitude near the stellar surface, indicating 
that the sensitivity to the rotation is the highest close to the surface. 
This reflects the fact that acoustic waves, the main constituents of 
high-order p modes, slow down as they come close to the surface, 
so that the mode periods, which are essentially determined by the 
residence time of the waves at each layer, are most influenced by 
the near-surface structure. Note that at the outer edge of the con¬ 
vective core of Kepler 25 and HAT-P-7, a small kernel discontinu¬ 
ity appears. This is caused by those of the chemical composition at 
the top of the convective core and is due to the fact that diffusion is 
neglected in our models. 

Figs Eb) and l^b) show the cumulative integrals 
fg°K„^i(r)dr for degrees Z = 1 and I — 2 , respectively. 
Note that within the observable range of n, which is approximately 
between 10 and 25, all of the kernels with I = 1 and I = 2 are 
nearly identical to each other. Still, there is clear dependence 
of the base of the convective envelope on the mass, such that 
more massive stars have shallower convective envelopes. This 
implies that the contribution of the radiative zone to the rotational 
splittings is greater in more massive stars than in less massive 
ones. Furthermore, the profiles in the radiative zone are almost 
linear functions of radius, between ~ 0.15i? and the outer edge of 
the zone. Since the gradient of the profiles indicates the sensitivity 
to each layer, each observed eigenmode probes the radiative zone 
nearly uniformly from the very deep inner layers to the outer edge. 


2.3 Rotation rates in the envelope and in the radiative zone 


A solar-like star (slightly more massive than the Sun) consists of 
an outer convective zone, an intermediate radiative zone and a core 
where nuclear fusion occurs. Figs [T] and demonstrate that low- 
degree modes are insensitive to the stellar core. In order to probe 
the rotation rat e within the (radiative) core, it is required to measure 
g modes (e.g. iGarcIa, Mathur & Ballol|[200i: iMathur et alj|2008l: 

lAppourchaux et al.l2010l). which are still elusi ve in solar-like stars, 

or mixed modes (e.g. Deheuvels et ahlEoidh . which are seen in 
evolved solar-like stars only. However, we stress that the rotational 
splittings of low-degree p modes in main-sequence solar-like stars 
are sensitive to the rest of the structure, both of the radiative zone 
and the convective envelope. 

In order to separate the contributions from the radiative and 
(outer) convective zones, it is assumed that the two zones rotate 
uniformly with different rates along the same axis. Then, it is pos¬ 
sible to express the rotational splitting, Sun^i, by a weighted aver¬ 
age of the rotation rate of the radiative zone, /rad, and that of the 
convective zone, /conv, as 


— trad frad “f Fconv /conv, 


(3) 


where Rad = K„,i{r) dr and 7conv = Kn,i{r) dr denote 

the integrals of the kernel in the radiative and convective zones, 
respectively, with Vc being the radius at the base of the convective 
zone. Note that 


Fconv “t” Rad — 1, (4) 

because each kernel is normalised to 1 . 

It is required to make two additional assumptions in order 
to make use of equation Firstly, the rotation rate in the con¬ 
vective envelope, /conv, is assumed to be equal to that at the 


surface, /surf. This is supported by the solar case, which shows 
30 per cent difference at most in the rotation rate in the latitu¬ 
dinal direction and m uch smaller in the radial direction in the 
convection zone (e.g. Sdhou_et ^ 1 19981: [Thompson et'ZI l2003l: 
lEff-Darwich & Korzennikl 20131) . Secondly. it is assumed that the 
rotational splittings remain nearly constant over the observed 
ranges of n and Z, so that the average rotational splitting, {Svnf, 
can be used as a representative value of the seismically measured 
internal rotation rate, /seis. This point has been justified in sec¬ 
tion | 2 j^ by demonstrating that the integrals of the kernels of the 
observed p modes differ little from each other. We can now use 
equation (0 to express the relation linking the surface rotation rate, 
/surf, the average rotation rate in the radiative zone, (/rad), and the 
average rotation rate of the stellar interior, /seis = (Si'n,i), as 


{frad) — /surf “t“ {Rad) (/seis /surf). 


(5) 


Note that () denotes the average values over the range of observed 
n and Z, with 10 < n < 25 and 1 ^ Z St 2. 

If we go back to equation 0 , it tells that /seis is reduced to a 
simple average of Rad and /conv, if {Rad) = (/conv). However, in 
stars with a thin convective envelope {i.e. the most massive solar¬ 
like stars), the sensitivity of the splitting to the rotation within the 
radiative zone may dominate, which means {Rad) (/conv), so 
that these stars are ideal candidates to probe a differential rate be¬ 
tween outer and inner layers. In fact, equation 0 is equivalent to 
/seis - /surf = {Rad) {{Rad) “ (/conv)), whose left-hand side is 
an observable. This equation shows that, for a given difference in 
the rotation rate between the radiative zone and the convective en¬ 
velope, the observed difference between Reis and /surf is larger (in 
magnitude) in stars with larger {Rad)- 

Applying the error propagation law to equation 0, we esti¬ 
mate that an interior with fas ter-than-surface rotation can be de¬ 
tected at a 1(7 detection level if (/seis — /surf)//surf ^ 25 per cent, 
which corresponds to (/rad)//surt sj 1.5. Similarly, detecting a 
slower-than-surface rotation requires (/seis — /surt)//surf R, —45 
per cent, or (/rad)//surf < 0.2. This is obtained by assuming that 
the relative uncertainty on /surf is of 30 per cent, based on the Sun. 
Here, the uncertainty account for the fact that /conv and /surf may 
not be identical, depending on the latitude at which the star is ob¬ 
served (latitudinal differential rotation^ We also set the uncertain¬ 
ties on /seis to 5 per cent and to 4 per cent on {Rad) (justified in 
section|4l(. 

In the following, equation 0 is used to measure the rotation 
rate in the radiative zone. 


3 DERIVING THE STELLAR ROTATION 

This section describes how the surface and internal rotation is de¬ 
rived from the observables. 


3.1 Surface rotation rate 

The surface rotation rate /surf can be estimated using either (1) ro¬ 
tational broadening of spectroscopic absorption lines, which gives 
V sin i — 2nR sin i, with the radius, R, and the inclina¬ 
tion angle, i, determined by theoretical models and asteroseismic 
analysis, respectively, or ( 2 ) the brightness modulation of period 


® Actually, this is a very conservative value as our formal uncertainties on 
/surf are in fact of Ri 10 per cent. 
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Figure 1. The figures show the rotational kernels (a,b) and their integral 
over the radius (c) for the Sun, Kepler-25 and HAT-P-7 for I = 1 modes. 
Plots have two scales to magnify the Kernel of near-surface layers. Only 
modes with radial orders n = 10 (a) and n = 25 (b) are shown. Inte¬ 
grals of kernels are almost identical to each other independently of the star. 
Vertical dotted lines indicate bases of the outer convective zones. Small dis¬ 
continuities around 0.1 in (a) and (b) are caused by those of the chemical 
composition at the top of the convective core. 


Figure 2. Same as Fig.^ but for 1 = 2 modes. 


(‘ 2 ') 

Prot = l//surf due to structures at the surface of the star (e.g. 
magnetic spots). Note that the rotational broadening is determined 
by contributions from a wide range of the area on the visible hemi¬ 
sphere, while the brightness modulation measures the rotation pe¬ 
riod at the latitude where spots evolve. For the Sun, it is well known 
that the spots migrate from « 30° llMaundeijll904h toward the 
equator during their 11-year cycle, so that v sini and Prot may 
not be sensitive to rotation of the same regions. This means that 
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the observed rotation frequency and may intrinsically 
differ from each other in case of significant latitudinal differential 
rotation. This point is further discussed in section|4] 

The spectros copic v sinf o f the analysed CoRoT stars has 
been measured by iBruntJ ( l2009h . while, except for HAT-P-7 and 
Kepler-25, t he studied Kepler s tars ha ve been observed s pectro- 
scopically bv iBruntt et alj ll20i^ . Since iBruntt et alJ j2012h do not 
provide uncertainties, we assumed here a typical precision of 10 per 
ce nt. Regarding HA T-P -7 and Kepler-25, w e use values determined 

bv |Paletal.U2008ll and iMarcv et all (l2014lt . respectively. 

( 2 ) 

The surface rotation frequencies of the two CoRoT 

stars, HD 49933 and HD 181420 are from iBenomar et all j2009lt 
an d iBarban et alj l l2009h . respectively, who have detected the low- 
frequency periodicity in the power spectrum. Note that alternative 
values were derived by using a stellar spot modelling technique in 
iMosser et alj ( l2009h . These are in agreement with the surface rate 
used in this study (see Table lA^ . For Kepler stars, they are de- 
termined from rota tion periods Prot, which have been measured by 
lOarcta et al.l ( 12014h based on the wavelet transform and on the anal¬ 
ysis of the autocorrelation function of the light curvfl When the 
comparison is possible, their derived surface rotation is consistent 
with the rotation from iKaroff et alj j2013h . Note that cannot 
be evaluated for HAT-P-7 and Kepler-25 because they do not show 
detectable surface activity. 


3.2 Internal rotation and stellar inclination 


Asteroseismology of solar-like stars, in the present case, is based on 
the photometry of the luminosity variation. Because of its stochas¬ 
tic nature, each solar-like mo de has a Lorentzian profile in the 
power spectrum ( lHarveviri985h . In that case, the stellar oscillations 
can be expressed as a sum of Lorentzians over n, I and m. 


= E E T 


H, 


n,l,m 




(6) 


Here, is given by equation ([T} , while and rn,i,m 

correspond to the mode height and width at half maximum, respec¬ 
tively. 

The turbulent convection, which is the origin of acoustic 
modes, does not prefer any particular direction to others, so that 
the rotationally split modes with the same I and n are expected to 
have a lmost the same amplitudes. Due to a geometrical projec tion 
effect dBallot. Garcta & Lambert|[2006l : lGizon & SolankJlTOoHl . in 
disk-integrated photometry, the height of the split modes depend on 
the stellar inclination angle i as Here Hn,i 

is the intrinsic height for the mode, while is the visibility 

of the m-components in the power spectrum. 




[I - |n^l)! 

{I -b |m|)! 




2 


(7) 


with being the associated Legendre function. It is therefore 
possible to determine i from Several authors have al- 

ready exploited that approach to measure the stellar inclination (e.g . 


Ballot et al.ll2008l: 

Benomar et al.ll2009l;lADt)ourchaux et al.ll2012l 

Chaolin et al.ll2012 

:lBenomar et al.ll2014al; Lund et alj|2014). 


In practice, the rotational splittings and the inclination an¬ 
gle are determined, together with the Lorentzian parameters, by 


Using KADA CS and PDC-MAP data, see lOarcta et at ] j2Qllh : 
iThompson et al.l i2Q13h . 


fitting the power spectrum of the light curve with the formula 
given by equation l|^. Therefore, the precision on the rotation is 
limited by the correlation between the rotational splitting and of 
the stellar inclination. Note that the splitting is assumed constant 
over the observed rang e of modes in this analysis. As shown by 
IBenomar et al.l ( 1201 4alt . the correlation problem is particularly im¬ 
portant in slow rotators, because the 2/4-1 m-components of 
the same order n and degree / are difficult to disentangle from 
each other. Contrary to heat driven pulsators, stochastically driven 
pulsators have an important intrinsic noise of stellar origin that 
renders mode pulsations hard to analyse and requires sophisti¬ 
cated statistical method s. Here, we use a Markov Chain Monte 
Carlo method (see e.g . Isenomar. Appourchaux & BaudinI 1 20091 : 

iHandberg & Campanta 1201 iL for applications in asteroseismol- 
ogy). This allows us to evaluate the probability density function 
for the rotational splitting and stellar inclination, as well as their 
correlations. Note that par t of s tudied stars have been ana lysed by 
lAppourchaux et al] ( l2012ll and lAppourchaux et al] ( l2014h . whose 
work focuses on mode frequencies, widths and heights. Stars of the 
current work in common with these studies are revisited here using 
the same data, but focusing on the extraction of the rotation and the 
stellar inclination angle. 


3.3 Centrifugal distortion 


In section im we have assumed that the internal rotation produces 
equal spacings in each multiplet of the power spectrum, as the 
first-order perturbation analysis due to the Coriolis force predicts. 
However, this cannot be fully justified if the second-order effect 
caused by the centrifugal force affects the pulsation modes signif¬ 
icantly. Assuming uniform rotation, we can estimate the effects of 
the Coriolis and centrifugal forces by the two dimensionless pa¬ 
rameters, ei = f /v and 62 = (//i^k)^, respectively, in which / 
and u are the rotational and pulsatio nal frequen cies, respectively, 
andz/K = {2n)~^(GMR~^Y^^ (e.g. lSaid[l98lh . 

These two parameters are related to the asymmetry of each 
multiplet by 



2 J L 2m J 


£2 

61 


( 8 ) 


for m ^ 1, where Vm represents the frequency of the component 
specified by the azimuthal order m. It is possible that the Coriolis 
and the centrifugal force are of the same order in some of the solar¬ 
like stars, which would induce an asymmetry of the multiplets in 
the power spectrum. Still, fitting such asymmetric structure with the 
symmetric formula should mainly introduce additional systematic 
errors in our estimates of the rotational splittings (and the stellar in¬ 
clination), rather than lead to completely wrong values of them. In 
fact, the frequency difference, {v,n — v-m)! (2m), which is equal 
to the rotational splitting, is not modified by the centrifugal distor¬ 
tion, but only changed by the third- and higher-order effects, which 
are negligibly small in all the stars. Moreover, if the asymmetry is 
significant, our analysis should suggest multiple different values of 
the rotational splittings. In the case of dipole modes, for example, 
we should find not only the correct value, — zzm=-i)/2, but 

also two other different spacings (at most), one between m = 1 and 
m = 0 components, and the other between m — 0 and m = — 1 
components. Since most of the probability density functions, which 
are shown in Fig. lAll do not have such multiple-peak structure, 
we presume that the asymmetry is generally smaller than (or com¬ 
parable to) the width of the main peak in the probability density 
functions. 
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3.4 stellar models 

In order to evaluate /conv and Irad in equation l[3]l, it is necessary 
to determine the thickness of the convective and radiative zones. In 
addition, an estimate of the radius of each star is needed to eval¬ 
uate the surface rotation rate based on the usini measurements. 
All of these quantities are obtained by computing stellar models 
that simultaneously match non-seismic observables (Teg, log g and 
[Fe/H]) and seismic observables (mode frequencies). The used 
constraints, their source and our results are summarised in Table 

lAll 

The best fitting models have been found using the ‘as- 
tero’ module of the Modules for Ex periments in Stella r Astro - 
physics (MESA) evolutionary code ( IPaxton et alj 1201 ll . l2013h . 
Stellar models have been calculated assuming a fixed mixing- 
length parameter aMLX = 2.0. We have used OPAL opacities from 
llglesias & Rogers! ( Il996h an d the solar compositio n (for the mix¬ 
ture of heavy elements) from lAsplund et al] i2009h . The initial hy¬ 
drogen abundance has been fixed to Xo = 0.7. No diffusion has 
been taken into account. Nuclear reactions have been set to include 
standard hydrogen and helium burning; the pp-chain and the CNO 
cycle in addition to the triple alpha reaction, which is the default 
setup of the MESA ‘astero’ module. We have utilised the NACRE 
compilation of nuclear reaction rates l lAngulo et alJll99^ with the 
update for ^“^Nfp, 7)^ ^0 and ^^C(a, 7)^®0 reactions iKunz et ^ 
l2002l;lFormicola et al.ll2004 . 

We have adopted the M. Schwarzschild treatment to define the 
boundary between the convective and radiative zones. Note that we 
do not take account of the convective overshooting, which we un¬ 
derstand modifies the age of the stars, the position of the Base of 
the Convective Zone (BCZ) and the stellar radius to some extent. 
In order to assess the importance of this effect, we have computed 
models of HD 49933 (index 21), HD 181420 (index 22), Kepler- 
25 (index 19), and HAT-P-7 (index 20) including the overshoot. 
These models have shown that by neglecting the overshoot, the 
BCZ and the radius are changed by only a few per cent, at a level 
that does not significantly influence the inferred internal rotation 
rate (/rad). However, it can affect the age (and hence the central 
hydrogen abundance) by 25 per cent typically. 

Eigenfrequencies have been calculated using ADIPLS 
dChristensen-Dals gaard 1200 Sl) . which assumes adiabaticity. We 
have applied s urface effect corrections to the frequencies, fo l lowin g 
the method of iKieldsen. Bedding & Christensen-DalsgaardI l l2008l) . 
The search for th e best fitting model ha s involved a simplex minimi¬ 
sation approach l lNelder & Meajl965tl using the criteria, where 
the mass M, metallicity [Fe/H] and age are the free parameters of 
the fit. 


4 RESULTS 

The sensitivity of the average rotational splitting to the rotation rate 
within the radiative zone is estimated to be (/rad) = 32-64 per 
cent (the full list is given in Table lAll in appendixj^. This range is 
consistent with what is reported in Figs[TJb) and[3b). Using the un¬ 
certainty in the radius, we estimate the error in (/rad) to be approx¬ 
imately 4 per cent. This is used to calculate the total uncertainty in 
the average rotation rate of the radiative zone. 

Fig. El shows the HR diagram of our samples (blue dots). 
The numbers identify the stars and correspond to those in the in¬ 
dex column of Tables lAll and IA2I Luminosities of these mod¬ 
els are determined so that the models fit both the non-seismic 



Figure 3. Best fitting models of the 22 analysed stars in the HR diagram. 
The numbers identify each of stars and con'espond to those in the index 
column of Tables Al and A2. These models have been computed with the 
MESA code so that they match both the non-seismic observables (Teff, 
logg and [Fe/H]) and the seismic obseiwables (mode frequencies) simul¬ 
taneously. Evolutionary tracks are for the solar' metallicity. 

and seismic observables, while the effective temperatures Ten are 
those from spectroscopy. The coloured lines are evolutionary tracks 
with the adopted physics and for the solar metallicity. Our sample 
contains predominantly F stars — some of which have been ex- 
tensively ana l ysed for activity variations and magnetic cycles by 
iMathur et alJ ll2014h —, with masses spanning from ~ 1.0 Mq to 
~ 1.6 Mq, an effective temperature 5990 K < Tefj < 6690 K and 
a metallicity similar to the Sun. Table IaTI provides more detailed 
information about the best fitting models of each star. 

4.1 Mode identification 

A fundamental step in asteroseismology is to fix the indices n 
(radial order), I (degree) and m (azimuthal order) of each ob¬ 
served mode. This is generally called mode identification. Eigen- 
modes of cool G type stars among our sample can easily be 
identified (at lea st for Z) by visual inspec t ion o f the so-called 
echelle diagram dOrec, Eossat & Pomerantzl Il983l) . Furthermore, 
the relative height of each low-degree mode in the power spec¬ 
trum (mode visibility Vi of degrees I, defined by H„^i = 
H„,i=o) is mostly due to a geometrical effect, and therefore 
is nearly independent of the fundamental p roperties of the stars 
dBallot. Barban & van’t Veer-Menneretl201lh . This property is also 
used to identify modes in the star’s power spectrum. 

However, in some cases, particularly for relatively hot stars, 
eigenmodes of the same parity (e.g. I = 0 and I = 2) are hard to 
disentangle from each other by visual inspection, mostly because 
the mode width at half maximum [see equation lU for its definition] 
is of the same order as the frequency spacing between the modes. In 
such a case, the most likely mode identification relies on a statisti¬ 
cal criteria such as the Bayes factor, between competitive solution^ 

® There ai'e actually only two possible solutions in the case of main se¬ 
quence solar-like stars. Readers are recommended to refer to the papers 
cited above for additional explanations. 
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lAppourchaux et a~ 20121 : ICorsaro & De RiddeJ 2014h . This prob¬ 
lem of mode identification is actually recurren t in F stars and was 
first e ncountered in the CoRoT star HD 49933 jAppourchaux et al] 
l2008h. In this paper, the mode identification of Kepler stars is based 
on I Appourchaux et alJ ( l2012h . which used the methods discussed 
above. As for HD 49933 (index 21) and HD 181420 (index 22), 
the identification reli ed on the Bayes factor. It is the same as in 
iBenomar et al.l (l2 009li fo r HD 49933, and corresponds to Scenario 
1 of Barban et alj l 2009h . for HD 181420. These identifications are 
also found as the most likely bv lWhite et alj l l2012h . using an inde¬ 
pendent approach. 

Interestingly, an incorrect mode identification of I will often 
lead to a disagreement between the spectroscopic v sin i and the 
seismic Useis sin i — 2nR /aeis sin i (see Table lA2l ). Therefore, 
stars with indications of important disagreement should be care¬ 
fully discussed. Here, we check that the mode identification is cor¬ 
rect for KIC 9139163 (index 10) which is also discussed in section 
14.61 We note that the modes of the same parity are disentangled and 
visible in the echelle diagram, with the expected mode visibilities. 
We therefore conclude that our current mode identification of I is 
correct for this star. 


As for KIC 3424541 (ind e x 3), the mode identification 
reported in lAppourchaux et aP ( 1201 2h leads to the seismic 
Uaeis siuI = 12.1 ± 2.3kms“^, which is incompatible with the 
spectroscopic v sini = 29.8 kmfrom iBruntt et alj (l2012h . In 
order to check the mode identification, we have used the latest Ke¬ 
pler data (Q5 to Q17), computed the Bayes factor of the two possi- 
bl e solutions and found that it is not in favour of the identification 
of lAppourchaux et H] ( l2012h . Note that due to its relatively fast ro¬ 
tation this star has a rotational splitting of the same order as the 
frequency spacing between modes of degree I = 2 and / = 0, 
while the stellar inclination is of « 90°. In such a configuration, 
the mode identification is not obvious because the azimuthal order 
m = —1 of the degree I = 1 can easily be falsely identified as 
an Z = 2 multiplet, while the m = -\-l could be falsely i denti- 
fied as an Z = 0. Note also that lAppourchaux et ahl ^20121) used 
a Bayesian framework, but with a shorter dataset. The new iden¬ 
tification, which has permuted I = 1 and the pairs of Z = 0, 2, 
leads to Vaeis sini = 26.89)'lQ g 4 kms~^. This new result is in 
agreement with the spectroscopic v sini. Additional information 
regarding the analysis of this star is available in appendix I bI 


4.2 Probability density functions 

Fig. a shows the joint-probability density functions of the stel¬ 
lar inclination and the rotational splitting, derived from asteroseis- 
mology for KIC 7206837 (index 8), KIC 9139163 (index 10) and 
KIC 9206432 (index 11). We also estimate the projected equato¬ 
rial radial velocity from this diagram by assuming uniform rota¬ 
tion; Vaeis siii i — 2 tiR faeis siii i. It also shows the surface ro¬ 
tation rates derived from two different kinds of observations as a 
function of the stellar inclination i; one from the lightcurve mod¬ 
ulation due to spots, and the other from the rotational broadening 
of spectroscopic absorption lines. KIC 7206837 is an example of 
good match between spectroscopic v sini, Prot, fseis and i, while 
KIC 9139163 and KIC 9206432 are examples of mismatch among 
the observables (candidates for significant radial differential rota¬ 
tion). The joint-probability density functions for all of the 22 stars 
are shown in Fig. lAll in appendix I aI 


4 


6 8 10 12 20 40 60 

Vggis sin i (km.s"’) Stellar inclination (degree) 

(a) KIC 7206837 


80 




6 8 10 12 14 15 20 25 30 35 40 

Vggis sin i (km.s'’) Stellar inclination (degree) 

(b) KIC 9139163 


Figure 4. The colour maps show the correlation between the stellar internal 
rotation rates measured by the rotational splitting and the stellar inclina¬ 
tion angles of three stars, on which superimposed are the surface rotation 
rates derived from the stellar spot (horizontal black dash-dotted lines) and 
the spectroscopic v sin i (curved dotted lines). The confidence intervals are 
shown in grey. In case (a), the surface and internal rotation rates as well as 
the stellar' inclination angle ai'e all consistent with each other, which indi¬ 
cates a near solid-body rotation. In case (b), although the indicator of the 
surface rotation is consistent with the seismically derived stellar inclination, 
it is in disagreement with the internal rotation rate. 


4.3 Surface rotation 

As the surface rotation rate depends on v sin i, the radius R 
and the stellar inclination i, it is useful to assess their respective 
contribution to the total uncertainty of We estimate the uncer¬ 
tainties of the radius, v sin i and the inclination i to be 5R/R ~ 1 
per cent, Sv sin i/v sini < 10 per cent and S sini/ sini ~ 16 per 
cent, respectively (see Table lA^ in appendix^. Thus, the uncer¬ 
tainty of the stellar inclination is the dominant error term in 
while the uncertainty in v sin i is the second most important term. 
Meanwhile, as for our estimates of stellar radii listed in Table lAll 
the results of 91 per cent of the stars are compatible within 2cr with 
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Figure 5. Comparison of the surface rotation rates estimated from (1) the 

rotational broadening of spectroscopic absorption lines, and (2) the 

( 2 ) 

lightcurve modulation due to spots, /g^^f Colours show the mass of stars, 
while the black dotted line shows a one-to-one relation. The indices for the 
stars coiTespond to those in the first columns of Tables [An and lA2l 

the estimates made bv iMetcalfe et an( l2014lfl This fraction drops 
to 72 per cent when comparing with stellar r adii derived from the 
scalin g relation that is calibrated by the Sun jKieldsen & Bedding 
Il995h. Still, compatib ility at 3cr is kept— see further discussions 
in iMosser et alj ( 1201 3ll —. Even if we consider a 3a uncertainty in 
the radius, the main error sources remain the stellar inclination and 

V sini. 

Fig.I^compares the two independent estimates of the surface 
rotation rates, and which have been derived from the 

spectroscopic v sin i and the lightcurve modulation due to surface 
structures, respectively. Although they are overall consistent with 
each other, some stars deviate from a one-to-one relation by ap¬ 
proximately a factor of two. This could be (1) due to systematic 
errors in the spectroscopic v sin i, (2) because multiple spot clus¬ 
ters appear at the same time at the surface of the stars, whereas the 
spot method is accurate only if one main structure exists, or (3) 
due to an intrinsic difference between the two estimates that origi¬ 
nates from a latitudinal differential rotation. When small, the spec¬ 
troscopic V sin-i can be poorly estimated because of the growing 
contribution of the microturbulence and macroturbulence relative 
to the actual rotational broadening of spectral lines. Microturbu¬ 
lence is of the ord er of one kilometre per second (e.g. lBruntj2009l: 
iDovle et alj|2014h . i.e. significantly smaller than our v sin i and is 
not expected to be a major source of bias. Ho wever, macroturbu¬ 
lence increases with the effective temperature dValenti & Fischej 
l2005h . and reaches 4-8kms“^ in the case of F stars. As macro¬ 
turbulence broadens the spectral lines, neglecting it overestimates 

V sin i. This could explain that Fig.|5]shows a group of stars, such 
as KIC 10162436 (index 13) and KIC 10454113 (index 15), that lie 
below a one-to-one relation, but cannot explain the population of 
stars with ;$> such as KIC 6508366 (index 5) and KIC 


® They modelled all Kepler stars of our sample except Kepler-25 and HAT- 
P-7, by taking account of more physical effects than our models, such as the 
overshooting and the diffusion. 



Figure 6. Ratio of the surface rotation rates deduced from two different 
kinds of observational data (one from the rotational broadening of spectral 
lines and the other from the lightcurve modulation due to spots) as a func¬ 
tion of the stellar inclination. The indices for the stars correspond to those 
in the first columns of Tables|AT]and|A2] The linear fit that includes all stars 
-case (a)- is shown in dai'k gray. The best fit is the solid line, uncertainties 
at 1(7 are the dotted lines. 

10355856 (index 14). However, we remark that these stars are seen 
from latitude near the equator (60° < i < 90°), where most stellar 
spots are expected to emerge. Therefore, the important discrepancy 
between and could be due to the presence of multiple 
spots. 

A latitudinal differential rotation may also lead to discrepan¬ 
cies between and To estimate this, we looked at the 

ratio /surf//surf ® function of the stellar inclination i (Fig.|^, 
which exhibits a trend. To assess its significance, we fitted a linear 
function /surf//surf ~ CL + bi, in three different situations, 

(a) with all observed stars. We obtain a = 0.26 ± 0.12 and 
& = 0.0015 ±0.008, 

(b) excluding the star which shows the most extreme difference 
(star 14). We obtain a = -0.085 ± 0.15 and b = 0.001 ± 0.003, 

(c) excluding the three most extreme stars (star 13,14,15). We 
obtain a = —0.18 ± 0.15 and b = 0.000 ± 0.003. 

While case (a) indicates a significant slope (see also Fig. |^, case 
(b) and (c) show that the fit result strongly dependents on a few ex¬ 
treme observations. We conclude that our measures do not strongly 
indicate evidences of a trend with i that could explain discrepan¬ 
cies described above. Further study on a larger sample, that better 
accounts for the macroturbulence is required before drawing firm 
conclusions on the latitudinal differential rotation. In any case, pos¬ 
sible biases in v sin i are expected to be less important than those 
in Prot, and in the following, we choose to for the indicator 
of the surface rotation. 

4.4 Comparison between internal and surface rotation 

Fig.|7] compares with the seismically obtained rotation rates 
/seis. For reference, the Sun is also shown. The solar values of 
the rotational spli tting for Z = 1 and I = 2 with 14 < n < 23 
were taken from iToutainI l l200lh . We obtained their average of 
{5v„i)q = 0.370 ± 0.007 upHz. The surface velocity is calcu- 
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Figure 7. Rotational splittings versus surface rotation rates measured us¬ 
ing the spectroscopic v sin i for sampled stars. The solar symbol near the 
lower left comer indicates the Sun. The numbers attached to the stars corre¬ 
spond to those in the first columns of Tables lAn and lA2l The ratios of these 
two quantities are mostly close to unity, implying that there is no signifi¬ 
cant radial differential rotation in most cases. The dashed triple-dotted line 
represents the expected positions of ZAMS models that assume local con¬ 
servation of angular momentum in the radiative layers and constant rotation 
rates in the convection zone. 


lated for the Sun-as-a-star, using th e solar vq sin iq = 1.6 ± 0.3 
kms~^ from IPavlenko et al.U 2012 h and adopting R 0 = 6.96x10® 
km. The ecliptic plane is inclined from the solar equatorial plane by 
approximately 7°, which means that sin iq can be different from 1 
by 0.8 per cent at most. In Fig.Q /sets are in good agree¬ 

ment, which already suggests that there is not a significant radial 
differential rotation in most cases. We thus confirm that nearly uni¬ 
form rotation, which has been well established in the case of the 
Sun, extends to the faster rotating main-sequence stars with simi¬ 
lar masses. However, we note that KIC 9139163 (index 10) depart 
significantly from the one-to-one relation, indicating that it seems 
to have faster-than-surface interior rotation. 

In order to discuss the implication of this result, evolution of 
the rotation profiles of 1.2 and 1.5 Mq stellar models with the 
solar composition has been calculated (Saio, private communica¬ 
tion) from the (pre-main-sequence) Hayashi phase to the Zero-Age 
Main-Sequence (ZAMS) stage, with the following assumptions: (1) 
rotation rates are slow enough to have negligible effects on the 
structure, which keeps spherical symmetry, ( 2 ) total angular mo¬ 
mentum is conserved (no addition to and no removal from the sur¬ 
face layers are considered), (3) the rotation profiles are functions 
of only radius, (4) angular momentum is locally conserved in the 
radiative zone and (5) it is redistributed instantaneously in the con¬ 
vection zone to establish a constant rotation rate (a limiting case of 
the efficient transport) with the total angular momentum of the zone 
fixed. Note that point (5) leads to uniform initial rotation profiles, 
because the models are fully convective in the Hayashi phase. 

As the models evolve towards the ZAMS, the central mass 
concentration gets higher, which results in spin-up of the cen¬ 
tral part, whereas the rotation rate in the envelope augments more 
weakly. Consequently, the ratio between the central and surface ro¬ 
tation rates reaches about 2.5 (2.7) for the 1.2 (1.5) Mq model. 
Note that this ratio is independent of the initial rotation rate itself. 


By substituting the calculated rotation profiles of the ZAMS models 
into equation 13, we have computed the average rotational split¬ 
tings normalised by the surface rotation rate over all modes with 
( = 1, 2 and n = 10-25. The resultant ratios are equal to 1.22 and 
1.24 for the 1.2 and 1.5 Mq models, respectively, which suggest 
that the ratio increases as a function of mass very mildly (at least in 
the range that we discuss in the present analysis). To take account 
of the smallest mass in our sample (~ 1.0 Mq), a common ratio 
of 1.2 is adopted for the ZAMS models between 1.0 and 1.6 Mq. 
Fig.|7] shows a dashed triple-dotted line with this value of slope. 
During the subsequent evolution of the models after ZAMS, the 
central mass condensation continues, leading to an increase of the 
ratio as the models age. The line thus indicates the lower limit of 
the ratio that main-sequence stars in the mass range under consid¬ 
eration would have, if angular momentum were locally conserved 
in the radiative layers and fully redistributed in the convection zone 
to realise a constant rotation rate. Therefore, the stars that are found 
below the line in the figure provide clear evidence that angular mo¬ 
mentum has been transported inside the radiative zone and/or be¬ 
tween the radiative and convective zones of these stars. With a la 
confidence interval, the number of such stars is equal to 13 (those 
with indices 2, 3, 4, 6 , 8 , 9, 12, 15, 17, 18, 19, 21 and 22). At 
the 2a leveQ, it still remains 10 stars (2, 3, 4, 6 , 9, 15, 17, 18, 
19, 22)0 Because the same physical process should operate in the 
same physical conditions, we presume that angular momentum is 
efficiently transported inside solar-like stars in general. Note that 
the magnetic braking of the surface layers, which is su pposed to be 
effective for low-mass stars in our sample (lKraftll967h . is expected 
to increase the ratio of the rotation rates between the interior and 
the surface. The stars that have ~ 1.3 Mq or below, and found on 
or above the line on the figure (those with indices 10, 13, 14 and 
16) could reflect a significant effect of this mechanism. 


4.5 Interpretation based on the simple two-layer model 

To quantify the level of radial differential rotation in a different 
way, we evaluated the ratio (/j^d )/Isurf based on equation (|5](- Fig. 
[^shows that if we pickup randomly a solar-like star of our sample, 
there is 68 per cent probability that it has a rotation rate in the 
radiative zone not greater by 54 per cent and not slower by 41 per 
cent than the surface. Individually, the radial differential rotation 
does not exceed a factor of two, except for KIC 9139163 (index 
10). Although it is in principle possible to obtain negative rotation 
rates, the observables do not provide us with enough information 
to distinguish the direction of rotation. Therefore, we exclude the 
possibility of opposite rotation rates (dashed area in Fig. 


4.6 Particular cases of KIC 9I39I63 and KIC 9206432 

In this section, we further discuss the case KIC 9139163 and KIC 
9206432, which were more thoroughly analysed than other stars. 


^ The calculation has been done using the probability density functions and 
thus account for the skewness of the distributions. 

® However, the values of t;sin 2 of stars with indices 4 and 18 are below 
5kms“^, so that the contribution of macroturbulence could possibly be 
more significant than the estimated uncertainties. 
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4.6.1 KIC 9139163 

Our analysis of KIC 9139163 leads to values of spectroscopic 
V sin i which are less than half as small as the seismic Useis sin i. 
The spectroscopy gives v sini = 4.00 kms“^, while the seismi- 
cally derived value is UseisSini = 10.15 ± 0.95 kms“^. It does 
not seem that these meaningful differences can be fully explained 
by only the known biases in the spectroscopic v sin i (see section 
14.31 . The surface rotation derived from spots still appears to indi¬ 
cate an important radial differential rotation, as — 2.3 

(instead of /^adZ/surf — 3-8). In general, the mode identification 
of oscillations in F stars is harder than in G stars, and one might 
think that this could also be a source of discrepancy. However, that 
latter possibility is ruled out for this star (see section l4~T] for further 
discussion). 

A careful review of the assumptions in our analysis revealed 
that the centrifugal distortion may explain the apparent discrepancy. 
In our analysis and as discussed in section [53] this was neglected 
so that the internal rotation is assumed to produce equal spacings 
in each multiple! of the power spectrum. However, by comparing 
the two parameters for a representative frequency of the star, it is 
found that they can be of the same order. We thus performed a new 
analysis that includes the centrifugal dist ortion, in similar fashio n 
as it has been done for Kepler-410A by IVan Evlen et al.l ( l2014h . 
We found that the centrifugal force effect is negligible and did not 
modify neither the splitting nor the stellar inclinatiorfl We con¬ 
clude that second order effects on the splitting are not the cause of 
the discrepancy between the surface rotation and the internal rota¬ 
tion. 

Mode visibilities Vi^i and Vi ^2 which were considered as 
free parameters may also modify the results as an important cor¬ 
relation with the splitting may exist. This is because in hot stars, 
modes of degree I — 0 and 1 — 2 (see section l4Tt , potentially 
overlap. However, a new analysis that hxes them to the expected 
visualities = 1.5 and Vj ^2 = 0.5 did not modify neither the 
splitting nor the inclination. 

4.6.2 KIC 9206432 

For this star, the mode identification is not evident from the echelle 
diagram and requires to rely on th e Bayes factor, which is in favour 
of the mode identification from lAppourchaux et all (l2012h . Fur¬ 
thermore, the s ame identifi c ation is obtained using the empirical 
approach from IWhite et"^ ( 1201 2h . Note also that the less likely 
mode identification gives i ~ 55° and /sets — 3 pHz, which leads 
to a great discrepancy between spectroscopic and seismic v sin i. 
We thus assume in the follow ing, that correct mode iden tification 
is the one already presented in lAppourchaux et alj | |2012|) . 

The inferred rotation for the most likely mode identification 
is very sensitive to the fit assumptions. Our analyses revealed that 
the stellar inclination is relatively low, so that the m 7 ^ 0 com¬ 
ponents of the degree ( = 1 do not have a significant power 
and do not provide an accurate constraint on the rotational split¬ 
ting. Furthermore, the modes of degree I — 0 and 1 = 2 over¬ 
lap significantly, so that these can hardly be disentangled. With¬ 
out tight constraints on the mode visibilities (i.e. when mode vis¬ 
ibilities are free parameters) we noticed that the fit is unable to 
distinguish the small separation do 2 from the rotational splitting 
{Si/ni). This leads to visibilities far greater than expected. We found 

We found Siy^i = 3.16lpg® and i = 23.612 s- 



Figure 8. Left. Average radial differential rotation between the radiative 
zone and the surface as a function of the stellar inclination. Grey lines indi¬ 
cate a differential rotation of a factor of two. The dashed area indicate nega¬ 
tive rotation rates. Right. Probability density function for — 1 

of the sampled stars. The median is shown in green, while orange lines in¬ 
dicate the confidence interval at Icr. 

Viti = 2.2 ±0.3 and Ilia = 1.3 ±0.2 instead of V)ii = 1.5 and 
Vi ^2 = 0-5. It also gives a biased rotational splitting which leads 
to an apparent disagreement between spectroscopic and seismic, as 
V sin i = 6.70 kms~^ while Vseia sini = 13.66llg4 kms“^. 

Fixing VjLi = 1.5 and Vi ^2 ~ 0-5, reduce the degen¬ 
eracy of solutions and gives a better agreement (Useis sin i = 
8 . 21 II km s“^). We also evaluated the influence of Z = 3 modes 
in our analysis (with all visibilities fixed and VjLg = 0.07) and 
found an even better agreement (Vseia sini = 7.77 ZIq gg kms“^). 

Thus in that particular case, it is necessary to provide strong 
constraints on the mode visibilities in order to disentangle the 
power associated to the Z = 0 from its neighbour 1 = 2 and conse¬ 
quently, to accurately measure the rotational splitting. Note that in 
the following, we provide the results of the fit with fixed visibilities 
and that includes Z = 3 because the calculation of the Bayes ratio 
indicates that these are statistically significant^ 


5 DISCUSSION 

In the following, we further discuss the implications of our results. 

5.1 Sample selection 

The ensemble of 22 stars analysed in this paper are selected by the 
following criteria: (1) they are not known to be members of close 
binary systems; (2) they must show solar-like oscillations; (3) they 
must have a surface rotation velocity of at least a few kilometres per 
second to have accurate v sin i; (4) they are preferably hot solar-like 
stars (F stars), so that they have a thin convective envelope, a condi¬ 
tion that improves the accuracy on the measured rotation rate of the 
radiative zone. While (1) is an important criteria for having reliable 
spectroscopy, our sample representativity for the ensemble of solar¬ 
like population has to be discussed for (2, 3, 4). Among the initial 

By comparing the model that includes 1 = 3 with the one that does not, 
we found a Bayes ratio of Ri 70, which corresponds to a probability of 98% 
in favour of the model with 1 = 3.. 
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Spectral type 

Figure 9. Comparison of the v sin i in our sampte (btu e diamonds) with 
V sin i from the cataiog of iGtebocki & Gnacinskil J2005h . Upper and tower 
tines indicate the Icr intervai of cataiog popuiation. The middte tine is the 
median. 


2000 stars selected for asteroseismic observation (in short-cadence 
mode) in the Kepler field, only 500 showed sufficie ntly low activ¬ 
ity to enable the detection of solar-like oscillations dChaplin et alj 
I 2 OIII) . This may indicate that we are biased towards low-activity 
stars. Furthermore, the least massive and cooler solar-like pulsators 
are not represented in our sample mostly because of their weak, 
hard to detect, pulsations. Yet, we hypothesise that because the 
same physical processes should operate in all solar-like pulsators, 
our results might be extrapolated to cooler or more active stars. One 
may also argue that our sample is not representative in terms of 
V sin i. However no evident bias can be seen in Fig.|3 which com¬ 
pares the spectroscopic v sin i of o ur sample (blue symbols) wit h 
a catalog of v sin i for 28179 stars dGlebocki & Gnacinskill200^ . 
In this catalog, we only selected stars of K9 to AO type for which 
uncertainties were provided (over 9000 stars). Thus we think that 
our conclusion that most stars slightly more massive than the Sun, 
rotate nearly uniformly is valid for most of solar-like stars. Further 
analysis, based on a larger and homogeneous sample of stars would 
help to confirm this. 

5.2 Main results and their implications 

We have demonstrated that, for most of the 22 main-sequence stars 
in our sample, which have masses between 1.0 and 1.6 Mq, the 
surface rotation rate and the average rotation rate of the nearly 
whole star (except for the central 15 per cent region in radius) are 
almost consistent with each other. The difference between the two 
rates is so small for 10 stars that the simple evolutionary models 
without transport of angular momentum inside the radiative zone 
and between the radiative and convective zones are clearly rejected 
(see Fig.|7}. Expecting that the physical mechanism is universal, it 
is inferred that efficient transport of angular momentum generally 
occurs in solar-like stars. Assuming uniform rotation in each of the 
convective envelope and the radiative zone, it has been found for 
21 out of the 22 stars that the rotation rates of the two zones do not 
differ by more than a factor of two. 

As gyrochronology suggests that younger stars rotate faster 
than older stars, we attempted to identify signs of decrease of rota¬ 
tion rate in the interior {fff) by using two age indicators. 


(a) the hydrogen central abundance from the models. 

(b) the seismic observed ratio between the small separation 
do 2 = Vn-i,i =2 — tzn.i^o and the frequency spacing between con¬ 
secutive modes of same degree Azz = 

As shown by Fig. llOl a). the models of the stars have various ages 
since the core-hydrogen abundance spans between ~ 0.01 and 
~ 0.54, leading to X^/Xq ^ 1 per cent at the lower limit and 
Xc/Xq ~ 78 per cent at the upper limitH No trend suggesting a 
change of rotation rate along evolution is visible. This indicates that 
an angular momentum transport mechanism is efficient from the 
early main-sequence stage, and possibly even since the pre-main- 
sequence phase. 

Along evolution, models indicate that dog/Azz decreases, 
which might be associated to a decrease of rotation rate in low-mass 
stars with a significant magnetic braking. Instead, in Fig. llOl b). the 
rotation appear to increase with do 2 /Xu. Actually, because this in¬ 
dicator also decreases as the mass increase, the seen trend is mostly 
due to a mass-rotation relation. 

Interestingly, only one star, KIC 9139163 (index 10) seems to 
have a significantly faster rotation rate in the interior than the sur¬ 
face. We note that this star is among the youngest and massive stars 
in our sample, which may indicate that the near solid body rotation 
has not been established in these stars yet. The study of a larger 
sample, with possibly younger stars, may help to better evaluate 
when the angular momentum transfer mechanism becomes effec¬ 
tive and to derive timescales of such a process. On the other hand, 
additional analysis of the current candidate for a large radial differ¬ 
ential rotation is also required. For example, new high resolution 
spect roscopic observatio ns would be useful to confirm measures 
from lSruntt et al.l ( 1201 2h and to assess the uncertainties in v sini. 
KIC 9139163 could be an interesting case for the study of the dif¬ 
ferential rotation because it shows a large contrast between interior 
and surface rotation. However we stress that due to the low stel¬ 
lar inclination (~ 20°, which implies that the power of |m| > 0 
relative to the m = 0 is only of ~ 15% for ( = 1 modes and 
~ 35% for ( = 2 modes), our result could be very sensitive to 
the noise realisation and therefore, on the noise background model . 
This issue has also been discussed bv ICorsaro & De Ridden ( l2014h 
for that star. Thus, further studies are certainly required to confirm 
the unusually fast interior rotation of KIC 9139163. 

Because the mass range of our sample between 1.0 and 
1.6 Mq is about t he same as that of six s ubgiants and young red 
giants analysed by IPeheuvels et all l l2014h . the comparison of the 
two studies allows us to examine the evolution of internal rotation 
of stars in this mass range during the post-main-sequence stage. 
Since our analysis has demonstrated that these stars rotate almost 
uniformly in the main-sequence stage, the contrast (of about factor 
60 at most) in the rotation rate between the core and the envelope 
detected in the subgiants and young red giants should be under¬ 
stood as the outcome of the competition between the two effects, 
the efficient angular momentum transport and the expansion of the 
envelope afte r the main-sequence ph ase. This picture was already 
suggested bv IPeheuvels et all l l2012h . with the assumption of uni¬ 
form rotation at the end of the main-sequence stage. The present 
analysis strongly supports their assumption. 


Xo = 0.7 is the initial hydrogen abundance (see section [3i4l . 
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5.3 Implication in the problem of misaligned exoplanet 
systems 

Some of exoplanet systems show significant misalignment between 
the rotation axis of the host star and the normal of the orbital plane. 
It is even claimed that some exoplanets orbit in the oppo site direc¬ 
tion to the rotation of the h ost star (e.g. Xue_et al] 2014) . In order 
to explain this phenomenon, iRogers. Lin & Lai] ( 2012h proposed a 
hypothesis that the near-surface layers of the host star rotate in the 
opposite direction to the rest of the star, being caused by the angu¬ 
lar momentum transport by gravity waves that are generated at the 
outer boundary of the convective core. Since our analysis simply 
assumes that all layers of the stars rotate in the same direction about 
a common axis, the sign change of the rotation rate cannot be mea¬ 
sured directly by our analysis. However, if such reversal occurs, it 
should be still detectable by the analysis in this paper as a large dif¬ 
ference between the surface rotation rate and the average rotation 
rate over the large range of the star (the rotational splitting). The 
fact that most of the stars in our sample, including those with known 
exoplanets (HAT-P-7 and Kepler-25), rotate almost uniformly (see 
Fig.|7j clearly contradicts such reversal. Thus, the present analy¬ 
sis doe s not support the hypothesis proposed bv iRogers. Lin & Laul 
) l2012h as a general picture of F type main-sequence stars. 


5.4 Characteristics of the method 


Our method of measuring the rotation rates of solar-like pulsators 
is based on (1) usage of asteroseismic analysis and (2) combination 
with spectroscopic observations. In the following, these two points 
are explained within the context of the internal rotation measure¬ 
ments in asteroseismology. 

The first point is justified by the fact that the rotational split¬ 
tings detected in the main-sequence solar-like pulsators are sensi¬ 
tive to the rotatio n not only in the convective envelope bu t also in 
the radiative zone. lLund. Miesch & Christensen-Dalsgaarj ) l2014al) 
examined the rotation kernels of high-order and low-degree modes 
of a solar-mass main-sequence stellar model, taking account of de¬ 
pendence on not only radius but also co-latitude. They wrote in sec¬ 
tion 2.1, “It is clear for all these modes that the kernels are mainly 
sampling the outer parts of the star”. But this statement should 
be accepted with great care for general solar-like stars. While the 
contribution fr om the radiative zone amounts to about 37 per cent 
in their model ( iLund. Miesch & Christensen-Dalsgaar j2014bh . the 
value can rise to more than 50 per cent for more massive solar-like 
stars, as demonstrated by FigsIHb) and[3b) of this paper. We un¬ 
derstand that this is mainly because the convection envelope is thin¬ 
ner (in fractional radius) for more massive stars. The non-negligible 
contribution of the radiative interior means that the rotational split¬ 
tings observed in solar-like stars also contain a considerable amount 
of information about the rotation in the radiative zone. 

The second point reflects that it is difficult to distinguish 
the rotation rates between the radiative and convective zones 
based only on asteroseismic analysis. This is because the ob¬ 
served eigenmodes of each star have very similar rotation ker¬ 
nels to each other, as shown in Figs [T] and and therefore 
all of the rotational splittings provide almost the same weighted 
average of the inte rnal rotation rates. This has been known 
for a long time (e.g. lAerts. Christensen-Dalsgaard & Kurtzll2010l) 
and confirmed by recent space-based observations, for which 
only the average rotational splitting within the range of ob¬ 
served modes is robustly measured (e.g.lAppourchaux et al.ll2008l : 
[Ballot. Barban & van’t Veer-Mennerelll201 ll : Chaplin et alj|201^ . 




doj/Av 

Figure 10. Internal rotation as a function two age indicators: (a) the core 
hydrogen abundance Xc derived from stellar' modelling and (b) the reduced 
small separation do 2 / Atz. In (a), the vertical dot line indicate the initial hy¬ 
drogen abundance. Models do not account for overshoot which can modify 
Xc by 25 per cent. Selected stars cover a large pai't of the main sequence, 
but no evident slow down of the interior is observed along evolution. The 
trend seen in (b) indicate a stellar mass-rotation relation. . 


Recentlv. lNielsen et al.l ( l20l4) were unsuccessful in detecting vari¬ 
ations of the rotational splittings in Kepler stars. Therefore, an in¬ 
dependent source of information is required to examine radial de¬ 
pendence of the rotation rate. For this purpose, we have exploited 
the surface rotation rate measured by the spectroscopic v sin i com¬ 
bined with the average rotational splitting, allowing us to evaluate 
the degree of radial differential rotation. 
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APPENDIX A: SUPPLEMENTARY CONTENTS FOR THE 
MAIN RESULTS 

Provided in this appendix are the following pieces of supplemen¬ 
tary information about the analysis in the main part: 

(i) summary of the best fitting models in Table lAll 

(ii) summary of the rotation rates at the surface and the radiative 
interior in Table lA2l 

(iii) joint-probability density functions in Fig. lAll 


Table Al. Parameters used to constrain models (left), those der ived from models (m idd le) and mass/radius deriv ed from the scal ing rel ations (right). Effective temperatures ai'e from IPinsonneault et alj i?Q12h and 0^ 

the other non-seismic constraints ([Fe/H] and log^') are from (i) lBruntt et all i2Ql3) . fii) lMarcv et al.|j20l3) . (iii) |Pal et alJilQQSl) or (iv)i runtj l l20Q9h . The uncertainty in (/rad) is about 4 per cent. 


Name 

Spec, type 

T'eff 

Constraints 

[Fe/H] 

logs 

i/L© 

Results from models 

M/M0 /?/Rq Xc 

BCZ 

(/rad) 

Results from seal. rel. 

M/Mq /?/R0 


(MK) 

(K) 


(cgs) 



(frac. rad.) 

(%) 



1 

KIC 1435467(‘l 

F8 

6433 ± 86 

-0.11 ±0.06 

4.09 ± 0.08 

4.38 ± 0.08 

1.294 ±0.039 

1.686 ±0.020 

0.295 

0.825 

49.6 

1.47 ±0.04 

1.75 ±0.02 ' 

2 

KIC 2837475(‘l 

F5 

6688 ± 57 

-0.02 ±0.06 

4.16 ±0.08 

5.28 ±0.03 

1.553 ±0.028 

1.691 ±0.020 

0.462 

0.902 

62.8 

1.51 ±0.03 

1.69 ±0.01: 

3 

KIC 342454 iW 

GO 

6475 ± 66 

0.01 ±0.06 

3.82 ±0.08 

10.50 ±0.20 

1.556 ±0.018 

2.551 ±0.026 

0.213 

0.875 

58.1 

1.80 ±0.05 

2.67 ±0.02 1 

4 

KIC6116048(‘) 

F9 

6072 ± 49 

-0.24 ±0.06 

4.28 ± 0.08 

1.98 ±0.02 

1.134 ±0.027 

1.274 ± 0.020 

0.003 

0.730 

38.8 

1.02 ±0.02 

1.23 ±0.01 ' 

5 

KIC 6508366(‘) 

F6 

6499 ± 46 

-0.08 ±0.06 

3.94 ± 0.08 

7.46 ± 0.08 

1.436 ±0.021 

2.142 ±0.020 

0.249 

0.864 

55.9 

1.61 ±0.03 

2.22 ±0.02 

6 

KIC 6679371(‘) 

F5 

6598 ± 59 

-0.13 ±0.06 

3.91 ±0.08 

8.47 ±0.07 

1.457 ±0.021 

2.167 ±0.024 

0.275 

0.909 

64.2 

1.70 ±0.04 

2.30 ±0.02 

7 

KIC 7103006(‘) 

F8 

6421± 51 

0.05 ±0.06 

4.01 ±0.08 

6.04 ± 0.08 

1.417 ±0.025 

1.941 ±0.020 

0.295 

0.853 

53.8 

1.57 ±0.03 

2.00 ±0.01 

8 

KIC 7206837(‘) 

F8 

6392 ± 59 

0.14 ±0.06 

4.17 ±0.08 

3.92 ± 0.06 

1.307 ±0.027 

1.562 ±0.018 

0.494 

0.842 

52.0 

1.38 ±0.03 

1.59 ±0.01 

9 

KIC913915l(‘) 

G0.5(‘^l 

6134 ± 48 

0.11 ±0.06 

4.38 ± 0.08 

1.88 ±0.03 

1.206 ±0.077 

1.170 ±0.008 

0.554 

0.780 

44.3 

1.18 ±0.02 

1.16 ±0.01 

10 

KIC 9139163(‘l 

F8 

6405 ± 44 

0.15 ±0.06 

4.18 ±0.08 

4.08 ± 0.03 

1.378 ±0.028 

1.561 ±0.018 

0.452 

0.860 

55.0 

1.49 ± 0.03 

1.60 ±0.01 

11 

KIC 9206432(‘) 

F8 

6494 ± 46 

0.23 ±0.06 

4.23 ± 0.08 

4.18 ±0.05 

1.500 ±0.050 

1.551 ±0.010 

0.539 

0.878 

58.3 

1.57 ±0.03 

1.59 ±0.01 

12 

KIC9812850(‘) 

F8 

6407 ± 47 

-0.16 ±0.06 

4.05 ± 0.08 

4.94 ± 0.07 

1.343 ±0.013 

1.809 ±0.010 

0.274 

0.822 

49.1 

1.40 ± 0.03 

1.83 ±0.01 

13 

KIC 10162436(‘) 

F8 

6346 ± 108 

-0.08 ±0.06 

3.95 ± 0.08 

4.30 ± 0.03 

1.259 ±0.015 

1.971 ±0.018 

0.092 

0.735 

38.6 

1.45 ±0.05 

2.05 ±0.02 

14 

KIC 10355856(‘) 

F5 

6558 ± 56 

-0.19 ±0.06 

4.08 ± 0.08 

5.00 ±0.05 

1.273 ±0.016 

1.721 ±0.018 

0.290 

0.862 

55.3 

1.43 ± 0.03 

1.79 ±0.01 

15 

KIC 10454113(^1 

F9 

6197 ±47 

-0.06 ±0.06 

4.31 ± 0.08 

2.47 ± 0.02 

1.279 ±0.024 

1.285 ±0.018 

0.535 

0.811 

47.9 

1.22 ±0.02 

1.27 ±0.01 

16 

KIC 11253226(^1 

F5 

6690 ± 56 

-0.08 ±0.06 

4.16 ±0.08 

4.70 ±0.07 

1.355 ±0.016 

1.603 ±0.023 

0.430 

0.893 

60.7 

1.69 ±0.06 

1.73 ±0.02 

17 

KIC 12009504('l 

F9 

6230 ± 51 

-0.09 ±0.06 

4.21 ±0.08 

3.01 ±0.03 

1.140 ±0.074 

1.386 ± 0.008 

0.247 

0.851 

52.5 

1.28 ±0.03 

1.45 ±0.01 

18 

KIC 12258514(^1 

G0.5('‘l 

5990 ± 85 

0.04 ± 0.06 

4.11 ±0.08 

2.55 ±0.01 

1.068 ±0.034 

1.543 ±0.035 

0.010 

0.661 

32.0 

1.47 ±0.04 

1.64 ±0.01 

19 

&p/er-25(“l 

F(b) 

6270 ± 79 

-0.04 ±0.03 

4.278 ± 0.03 

2.64 ± 0.06 

1.285 ±0.036 

1.348 ±0.011 

0.500 

0.820 

49.0 

1.32 ±0.03 

1.36 ±0.01 

20 

HAT-P-7(“') 

F6 

6350 ± 80 

0.26 ±0.08 

4.070 ± 0.08 

5.84 ± 0.05 

1.540 ±0.030 

2.000 ± 0.009 

0.270 

0.870 

57.4 

1.48 ± 0.03 

1.97 ±0.01 

21 

HD 49933(‘''l 

F8 

6570 ± 60 

-0.44 ± 0.03 

4.28 ± 0.06 

3.76 ± 0.08 

1.287 ±0.009 

1.463 ± 0.006 

0.458 

0.880 

58.3 

1.34 ±0.02 

1.49 ±0.01 

22 

HD 181420(‘''l 

F6 

6580 ± 105 

0.00 ±0.06 

4.26 ± 0.06 

4.49 ±0.13 

1.335 ±0.010 

1.622 ±0.007 

0.400 

0.877 

57.8 

1.61 ±0.05 

1.73 ±0.02 


(a) Assumed as GO in Fig.[3 (b) Assumed as F5 in Fig. [9] 
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Table A2. Rotation rates at the surface of the stars and in the radiative zone. The penultimate column indicates the length of the data set while the last column gives the first and last used Kepler data quarters. 


index 

Name 

fseis 

i 

Wseis sin 2 

V sin 2 

±rot 

HI) 

^ surf 

A2) 

^ surf 

(/£d> 

(/i'd) 

Obs. duration 

Quarters 



(pHz) 

n 

(kms“l) 

(kms“l) 

(d) 

(oHz) 

( 9 Hz) 

(oHz) 

(oHz) 

(d) 


1 

KIC 1435467 

1 '70 + 0.23 
-*-•' "^- 0.20 

62 . 9 + 1 = 5 “ 

ii.zy _0 66 

10.00 ± 1.00 

6.68 ±0.89 

1 cq+ 0.22 

Y 9 + 0-23 
-*-•' '^-0.23 

l.yD-0.33 

1 70"*"0-50 

1- '0-0,47 

738 

Q5-Q12 

2 

KIC 2837475 

0 -|q+0.08 

o.iy_0.09 

75.4011“ 

22.7911“ 

23.50 ±2.35 

3.68 ±0.36 

q tio+0.34 

O.ZO_Q 33 

9 14+0-31 

■J-1^-0.31 

9 1 7+0.23 
'^•^‘-0.23 

q 9 c:+0.23 
O.ZO_Q 23 

738 

Q5-Q12 

3 

KIC 3424541 


84.9011“ 

26 . 8811 .®^ 

29.80 ±2.98 
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(a) Using spot modelling, ^losse^^L (200^ provide alternative surface rotation rate for the CoRoT stars. For HD 49933, they found 3.45 di 0.05 days and 2.25 ± 0.03 days for HD 181420. 
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Figure Al. Joint-probability density functions for all of the analysed stars. Horizontal dash-dotted lines indicate the surface rotation rates from lightcurve 
modulation due to spots (when available), while dashed curves are the rotation rates derived using the spectroscopic v sin i. Unless specified, all spectrum fit 
include modes of degree Z = 0,1, 2. 
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(7) KIC 7103006 



(8) KIC 7206837 



(9) KIC 9139151 






(10) KIC 9139163 



(11) KIC 9206432 (with / = 0,1, 2, 3 and fixed visibilities) 



(12) KIC 9812850 


Figure A1 - continued 
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(13) KIC 10162436 



(14) KIC 10355856 



(15) KIC 10454113 
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(17) KIC 12009504 
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Figure Al - continued 
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Figure A1 - continued 
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APPENDIX B: SUPPLEMENTARY CONTENTS EOR KIC 
3424541 

The best fit of the power spectrum for the two possible mode identi¬ 
fication of KIC 3424541 is shown in Fig. lBlI and Fig. lB2l The most 
probable mode identification was assessed by using the Bayes ratio 
0(Mi, M 2 ) between the two possible competing models Mi and 
M 2 (e.g. [Benomar. Appourchaiux & Baudinll200^ . 0{Mi, M 2 ) is 
the equivalent of the the likelihood ratio used in frequentist ap¬ 
proaches. However, contrary to the likelihood ratio, the Bayes ratio 
can be interpreted in terms of probability. For example a Bayes ra¬ 
tio of 0{Mi, M2)=3 in favour of Mi indicates that the model Mi 
has 75% of probability to be correct. It is often consider ed that a 
mode l is strongly supported when the ratio exceeds 100 jjeffrevsl 
Il96lh . In the case of KIC 3424541, the mode identification shown 
in Fig IBlI is very strongly supported because the logarithm of the 
Bayes factor is 7.5 x 10® in favour of this identification. 

The list of measured pulsation frequencies for the most likely 
mode identification of KIC 3424541 is given in Table lBlI 


APPENDIX C: SUPPLEMENTARY CONTENTS EOR KIC 
9206432 

Our thorough analysis of KIC 9206432 enabled us to identify the 
likely cause for the apparent discrepancy between the surface and 
the interior rotation. In Fig ICll (top), we show the results for the 
stellar inclination and the rotational splitting when the power spec¬ 
trum is fitted by considering the modes visibilities Vj^i and 
as free parameters. In that case, there is an apparent mismatch be¬ 
tween surface rotation and rotational splitting. When the visibilities 
are fixed to their expected values, the agreement becomes better 
(middle and bottom). 


Table Bl. Mode frequencies for the most likely mode identification of KIC 
3424541 presented in this paper. 


degree 1 

index 

frequency (|i.Hz) 

0 

1 

496.77 ± 0.58 

0 

2 

535.47 ±0.58 

0 

3 

575.46 ±0.34 

0 

4 

615.77 ±0.35 

0 

5 

658.44 ± 0.27 

0 

6 

699.47 ±0.28 

0 

7 

740.87 ± 0.32 

0 

8 

781.62 ±0.30 

0 

9 

822.33 ± 0.33 

0 

10 

864.54 ± 0.41 

0 

11 

905.98 ±0.52 

0 

12 

947.90 ±0.57 

0 

13 

991.82 ±0.67 

0 

14 

1035.99 ± 1.97 

1 

1 

474.25 ± 0.61 

1 

2 

514.14 ±0.38 

1 

3 

553.20 ± 0.30 

1 

4 

593.10 ±0.25 

1 

5 

634.54 ± 0.23 

1 

6 

676.47 ± 0.22 

1 

7 

718.29 ±0.23 

1 

8 

759.57 ±0.24 

1 

9 

799.76 ± 0.30 

1 

10 

841.50 ±0.36 

1 

11 

883.52 ±0.41 

1 

12 

925.42 ±0.56 

1 

13 

967.46 ± 0.55 

1 

14 

1009.04 ±1.16 

2 

1 

493.35 ± 1.67 

2 

2 

533.17 ±0.90 

2 

3 

573.26 ± 0.64 

2 

4 

614.03 ±0.53 

2 

5 

655.45 ±0.52 

2 

6 

697.51 ± 0.44 

2 

7 

738.99 ±0.51 

2 

8 

781.14 ±0.50 

2 

9 

822.45 ±0.55 

2 

10 

864.26 ± 0.63 

2 

11 

905.82 ± 0.76 

2 

12 

946.56 ± 1.09 

2 

13 

987.69 ± 1.46 

2 

14 

1028.68 ±2.59 
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Figure Bl. Fit of the power spectrum of KIC 3424541 with the most likely mode identification (opposite as in lAppourchaux et al.U20l3n . obtained from the 
Kepler data analysis using Quarters 5 to 17. In this identification, the first excess of power on the left of the figure is due to a pair of degree / = 2, 0, and is 
followed by alternating 1 = 1 and Z = 2, 0 modes. 
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Figure B2. Fit of the power spectrum of KIC 3424541 with the less likely mode identification (same as in I Appourchaux et alj j2012h t. obtained from the 
Kepler data analysis using Quarters 5 to 17. In this identification, the first excess of power on the left of the figure corresponds to a degree Z = 1, and is 
followed by alternating Z = 2,0 and Z = 1 modes. 
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Figure Cl. KIC 9206432. Comparison of the results obtained for the stellar inclination and rotational splitting by including; (top) modes of degree I = 0, 1, 2 
mode visibilities as free parameters; (middle) modes of degree / = 0,1, 2 with fixed visibilities; (bottom) modes of degree / = 0,1, 2, 3 with fixed visibilities. 














































































